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EFFECTS  OF  RADIATION  DAMPING  ON  BEAM  QUALITY 
IN  THE  INVERSE  FREE  ELECTRON  LASER  ACCELERATOR 


Introduction 


Electron  beam  quality  as  measured  by  the  transverse  emittance  is  usually 


determined  by  the  gun  and  propagation  configurations  in  accelerators.  Under 


idealized  conditions,  the  transverse  normalized  beam  emittance  remains  a 


constant  of  motion  as  the  beam  propagates  through  the  accelerator.  Therefore, 


to  improve  the  quality  of  the  beam,  it  is  necessary  to  decrease  the  beam 


emittance  at  the  injection  point.  However,  since  the  normalized  beam 


emittance  is  essentially  the  transverse  area  in  phase  space  for  the  collection 


of  beam  particles,  one  can  in  principle  reduce  the  emittance  if  a  dissipative 


mechanism  is  introduced.  A  natural  candidate  for  such  a  dissipation  mechanism 


is  the  induced  synchrotron  radiation  damping  due  to  the  transverse  motion  of 


the  particles  in  an  external  periodic  transverse  magnetic  field.  It  is  this 


mechanism  that  we  will  focus  on  when  the  external  magnetic  field  is  chosen  to 


be  a  helical  wiggler  field.  Since  this  radiation  damping  effect  is  small  at 


low  energies,  it  is  in  the  context  of  the  recently  proposed  high  energy  IFEL 


accelerators1  11  that  we  will  concentrate  in  this  paper. 


We  begin  by  obtaining  the  electron  orbits  in  an  IFEL  accelerator. 


A  fully  relativistic  formulation  of  the  equations  of  motion  which  include 


radiation  damping  force  is  considered.  The  damping  coefficients  are  obtained 


from  the  transverse  dynamics  of  the  particles  while  the  axial  dynamics 


describes  the  acceleration  of  the  particles.  In  the  second  section,  a 
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relativistic  envelope  equation  for  the  average  radius  of  the  electron  beam  is 


derived,  assuming  continuous  emission  of  the  synchrotron  radiation.  It  is 
apparent  from  this  envelope  equation  that  the  normalized  transverse  emittance 
decays  exponentially  ao  a  rate  given  by  the  radiation  damping  coefficient. 

The  envelope  equation  is  solved  using  a  WKB  method  in  the  third  section  and 
the  spatial  evolution  of  the  beam  radius  is  obtained.  Quantum  excitation  sets 
a  minimum  value  on  the  normalized  transverse  emittance  in  an  IFEL  accelerator 
and  it  is  derived  in  the  fourth  section.  Strong  focusing  is  found  to  be 
necessary  to  reduce  such  minimum  to  an  acceptable  value.  An  example  is  given 
in  the  last  section  for  a  set  of  proposed  IFEL  accelerator  parameters.  It  is 
found  that  radiation  damping  does  reduce  the  emittance  of  the  accelerated 
electron  beam  while  resulting  in  an  insignificant  loss  in  particle  energy. 

I .  Single  Particle  Dynamics 

We  shall  consider  the  motion  of  an  electron  under  the  influence  of  a 

helical  wiggler  field  and  a  circularly  polarized  electromagnetic  wave  with  the 

inclusion  of  the  radiation  reaction  force.  The  fully  relativistic  equation  of 
1  U 

motion  is 


d£  v  x  B 

dt  =  +  +  l 


where 


.£  Y2-  !  dlEl  2  d£  2  d  d£  , 

TR  ~~2~2  ^  Y  ^~dt_)  ”  Yvdt^  +  dt(Ydt^  ’ 

m  c 
o 


0  7  0 

is  the  radiation  damping  force,  tr  =  2je|  / 3m0c  ,  and  Y  =  1  +  |£|  /mQc  • 

The  radiation  field  is  given  by  its  vector  potential  A  =A  (cos$e  -  sin<j>e  ), 

Li  l*  x  y 

where  <J>  =  kz~wt.  We  shall  assume  z-dependence  for  both  the  magnitude  and 


period  of  the  wiggler  field.  The  vector  potential  A  for  the  helical  wiggler 

— w 

field  is  given  by  A  =A  [cosQe  +sin9e  j  where  A  =  A  (z)  and  9=1  k  (z')dz'. 
&  — ww  x  y  ww  ■'  w 

J  n 


The  requirement  that  the  wiggler  field  satisfies  both  V  •  B  =0  and 

— w 

7  x  Bw  =  0  introduces  transverse  variation  as  well  as  a  nonzero  z-component  of 
the  magnetic  field.  J 

Since  we  shall  be  primarily  interested  in  laser  driven  accelerators,  the 

p 

normalized  wiggler  field  strength  aw  =  |e|Aw/m0c  is  assumed  to  be  much 

O 

greater  than  the  corresponding  quantity  a^  =  |e|AL/m0c  for  the  radiation, 
i.e.,  a  >>  a.  .  It  can  then  be  shown  that  the  major  contribution  to  the 

W  L 

radiation  damping  is  from  the  wiggler  field. 

We  shall  first  look  at  the  radiation  damping  term  in  Eq.  (1).  By 

neglecting  the  transverse  dependence  of  the  wiggler  field  for  a  beam  that  is 

confined  sufficiently  close  to  the  axis,  we  have  the  immediate  consequence 

that  the  canonical  momenta  in  the  x  and  y  directions  are  constants  of  motion 

and  may  be  chosen  to  be  zero.  The  mechanical  momenta  are  then  given  by 

p  =  A  *e  ,  p  =  A  *e  ,  where  A_  =  A  +  A,  .  Also,  in  the  zeroth 
x  c  -T  x  y  c  -T  y  -T  -w  -L 

order  approximation,  the  total  relativistic  energy  is  conserved  which  leads 
to  Y  =  0  and  pz  =  0.  Therefore,  the  only  significant  term  remaining  in  the 
radiation  reaction  force  is 


R  d  £  £  d£  2 

-  =  TRYtTT2  2  2  (dt)  -1  * 

dt  me 

o 


Neglecting  terms  of  order  a. /a  <<  1,  the  components  of  the  radiation  reaction 

L  W 

R  R  R 

force  are  F  =  -  v,cp  ,  F  =  -  v.cp  ,  F  »  -  v.cp  ,  where 
x  1  Fx’  y  1  V  z  |  z 


v  =  iDYk2c(a2  *  1 ) , 
_[  R  w  w 


2  2 

v.  »  t_Yk  ca 


are  respectively  the  spatial  decay  coefficients  due  to  radiation  damping  in 

2 

the  transverse  and  axial  directions.  Note  that  v,  =  v.  for  a  >>  1  which  is 

1  I  w 

the  case  in  the  IFEL  accelerator. 

The  most  significant  feature  of  the  transverse  motions  of  the  electrons 
is  the  betatron  oscillation  caused  by  either  the  inhomogeneity  of  the  wiggler 
field  in  the  transverse  plane  or  other  focusing  mechanisms.  It  can  be  shown 
that,  for  small  oscillations  about  the  axis  of  the  wiggler  field,  the 
transverse  equations  of  motion  are, 

d2 x  v2  ,  Y '  .  dx  ... 

71  *  V  ‘  '<  7  '  V  Tz  ■  (3a) 

dz  x 

7§  *  kb*  ■  -<  7*  V  &  • 

where  d/dt  =  v  3/3z,  v  *  c,  '  *  3/3z  have  been  used,  and  Kn  is  the  wave 
number  of  the  longitudinal  betatron  oscillation.  For  betatron  oscillations 
that  are  originated  from  the  v  x  B  force  due  to  the  nonzero  magnetic  field  in 
the  z-direction  of  the  realizable  wiggler  field, ^  Ku  =  a  k  /(/2Y)  . 

D  W  W 

The  axial  motion  of  the  electron  is  governed  by 


-  (v  x  B)  -  v  pz  , 

c  z  * 


where 


~  |  e  |  v  •  E  (vj_pj_2  +  v  p2  ) 


2  3 
m  c  t 
o 


It  is  straightforward  to  show  that  the  axial  electron  acceleration  is 


dv  3a  2a  a,  k  c  2k  3v|a,  a 

z  cw  wLw.  w  -LLw 

—  =  -  — - * - J  Sln  * - v.vz  +  - 2  C0S  *  ’ 

2Y  3z  Y  k  1  Y 


where  to  =  e  +  <t  =  f  [k  +  k  (z  )  -  w/v  (z')]dz  +  41  is  the  phase  between  the 

J  w  z  o 

o 

electrons  and  the  ponderomotive  wave  generated  by  the  beating  between  the 
radiation  and  wiggler  fields,  and  is  the  initial  phase  at  the  entrance  of 
the  interaction  region.  Equation  (5)  can  be  transformed  into  the  following 
pendulum  equation 


dj(/ 

dz2 


dk 

_ V 

dz 


,  8a  2  a  a.  kk 
k  w  w  L  w 

+  - - -  sm  ip  - 

2Y2  3Z  Y2 


2v,k 

—L^.  + 

c 


3vlaLawk 


COS  Ip 


(6) 


CY 


The  rate  of  change  of  relativistic  energy  may  be  obtained  from  Eq.  (4)  and  is 


d Y 
dz 


a,  a  k 
L  w 


vj_a,  a 


sin  ip  - 


V 


-(£ - 2)  cos  ip 

w 


ll(a2  ♦ 

Y  w 


aL? 


(7) 


Equations  (3),  (6)  and  (7)  will  be  the  basic  equations  we  shall  use  in 

studying  the  effects  on  beam  quality  due  to  radiation  damping.  The  terms 

containing  cos  ip  in  Eqs.  (5),  (6)  and  (7)  as  well  as  the  last  term  in  Eq.  (7) 

2  2  2 

may  be  neglected  when  the  conditions  a  >>  a.  ,  a  >>  1 ,  k  >>  k  ,  and 

W  L  W  W 

2 

Y  >>  1  are  satisfied.  These  conditions  are  easily  achieved  in  high  energy 
IFEL  accelerators. 


II .  Derivation  of  Envelope  Equation  with  Radiation  Damping 

The  single  particle  equations  of  motion  that  we  have  developed  in  the 
last  section  will  enable  us  to  study  the  macroscopic  behavior  of  the  beam. 
This  is  accomplished  by  considering  the  evolution  of  various  averaged 

10  1-3 

quantities  associated  with  the  single  particle  variables.  We  begin  by 

multiplying  Eq.  (3a)  by  x'  and  x,  and  Eq.  (3b)  by  y*  and  y,  where  ' 
denotes  9/8z.  Combining  the  resulting  equations  yields  the  following  set  of 


equations 


2  dz 


2  2 
dz 


.  2  v2  2 

VV 


V  d  2 
—  —  r 
2  dz 


=  “U*.  , 


where  r2  =  x2  +  y2,  g2  =  x'2  +  y'2,  p  =  Y'/Y  +  v^,  and  £  =  (x'y  -  y'x)  is  the 

2 

normalized  angular  momentum.  We  eliminated  8^  by  substituting  Eq.  (8b)  into 
Eq.  (8a).  By  taking  transverse  ensemble  averages  over  beam  particles  in  Eq. 

p  p  p 

(8),  and  denoting  the  ensemble  average  of  r  by  a  =  <r  >,  we  obtain  an 
equation  which  governs  the  evolution  of  the  root-mean-square  radius  of  the 
electron  beam, 

2d  2  d2  2  ^  _  v2  2  d  ,p  d  2.  1  d3a2  ^  d  tv2  2. 

p  — -  a  +  p  — -  a  +  2pK  a  +  — (x  x-a  )  +  —  — --  +  — (K  a  ) 


,  a  +  a  +  2pkb  a 

dz 

+  K2  ^-a2  =  0  . 

B  dz 


dz  2  dz 


2.3  dz'  B 
dz 


It  is  easy  to  show  that  the  integration  factor  for  Eq.  (9)  is  g  a  where 

2  2  fZ 

g  =  y  exp(2j  v  dz')  .  Equation  (9)  can  now  be  put  into  the  form 
o 

2  o  o  jj  2 

d/dz[  g  (aJa''  +  paJa'  +  a  ) ]  =  0  ,  and  can  be  integrated  to  give 

D 

2  ^  3  4  2  2  2 

g  [a  a''  +  pa  a'  +  a  K  ]  =  H  ,  where  H  is  a  constant  of  motion  associated 

D 

with  the  beam.  It  can  be  shown  that,  using  the  following  representation  for 

1  2 

the  particles'  normalized  transverse  velocities, 


a'  "  Lr 
_  r  er  *  -  e9  *  SB  , 

A 


where  68^  is  the  normalized  transverse  velocity  spread,  and  L 

p 

Eq.  (8c),  the  constant  H  is  given  by 


=  <£>  from 


where  Y(0)  =  Y(z=0)  and  L(0)  =  L(z=0).  We  may  therefore  define  the  squared 
normalized  beam  emittance1 2 ’ 1  ^  as  e^(z)  =  Y2a2<  |6£^|2>  and  arrive  at  the 


following  envelope  equation 


2  ,  „  [e2(z)  +  Y2L2(z)] 

d  3  /  1  uY  x  d<i  2  n  ~ 

— o  +  (  v  T-  +  V.)  —  +  KDa  -  - r— r -  =  0 

.2  Y  dz  1  dz  B  2  3 


The  spatial  dependence  of  the  normalized  emittance  and  average  angular 


momentum  are  given  respectively  by 


e  (z)  =  e  (O)exp(-f  v,dz')  , 
n  n  J 


(11a) 


L(z)  =  (Y(0)/Y)  L(0)exp(~j  v  dz')  , 

o  -L 


(11b) 


where  e  (0)  =  cn(z=0).  Equation  (10)  together  with  Eq.  (11a,b)  constitute  the 
beam  envelop*,  equation  with  radiation  damping  terms  included. 

One  can  see  that  when  =  0,  Eq.  (11a)  shows  that  en  remains  constant 
and  Eq.  (10)  reduces  to  the  usual  relativistic  beam  envelope  equation  where 


e  is  the  familiar  normalized  beam  emittance. 
n 


12,16 


Therefore,  in  the  presence 


of  radiation  damping,  the  root-mean-square  beam  radius  is  still  described  by 
an  envelope  equation  but  the  normalized  beam  emittance  is  no  longer  constant 
but  decays  exponentially  according  to  Eq.  (11a).  However,  the  decay  of  the 
normalized  beam  emittance  will  eventually  be  limited  by  quantum  excitation  due 


to  the  discrete  nature  of  the  synchrotron  radiation.  It  is  shown  in  a  later 


section  that  when  an  equilibrium  is  reached  between  these  two  competing 


processes,  the  minimum  normalized  emittance  achievable  through  radiation 


damping  in  the  IFEL  accelerator  is  given  by  ( e  )  .  *3fiaJk  /(/2m  clO  . 

n  min  w  w  o  B 

In  the  presence  of  radiation  damping,  the  average  angular  momentum  also 
decays  exponentially  as  given  by  Eq.  (11b).  However,  one  may  choose  L(0)  =  0 
for  beam  generation  schemes  that  do  not  impart  an  average  angular  momentum  to 
the  electron  beam,  i.e.,  zero  magnetic  field  at  the  cathode.  We  shall  assume 
that  this  is  the  case  in  our  study  of  beam  quality.  We  shall  also  not 
distinguish  between  v  and  V|  ,  and  will  denote  both  by  v. 


Ill .  Evolution  of  Beam  Radius 

The  equation  for  the  root-mean-square  radius  a  in  Eq.  (10)  is  nonlinear. 

p 

It  is  found,  however,  that  the  mean  square  radius  a  satisfies  Eq.  (9),  which 
is  a  linear  differential  equation.  For  beam  focusing  provided  by  the  wiggler, 
Eq.  (9)  may  be  solved  exactly  for  untapered  wiggler  fields  when  Y '  =  0. 

If  Y’  *  0  or  when  the  tapering  is  known,  it  can  be  solved  using  a  WKB  method 
if  we  assume  the  coefficients  are  slowly  varying.  Equation  (9)  can  be 
simplified  in  certain  limits  of  accelerator  designs  to  facilitate  analytical 
study.  It  can  be  shown  that,  r'/Y  <<  KD  and  v  i  KD ,  which  allow  us  to  arrive 

D  D 

at  the  following  approximate  equation 


?uS' 


“kbs' 


♦  +2(K?)']S 


B 


B 


(12) 


p 

where  3  =  a  . 

In  order  to  obtain  net  acceleration  of  the  electrons  trapped  in  the 
penderomotive  potential,  the  wiggler  field  must  be  spatially  tapered.  In  such 
a  case,  the  envelope  equation,  Eq.  (12),  is  a  linear  differential  equation 
with  spatially  dependent  coefficients.  We  solved  it  by  using  the  WKB-method 
which  assumes  these  coeffients  to  be  slowly-varying  functions  of  longitudinal 


distance.  By  assuming  both  K'/K  and  p  <<  KD ,  the  general  solution  to 

Do  D 

Eq.  (12)  is  found  to  be 


-M  KB(0) 

S  =  e  — ■ 7  r  [  A  +  B  cos  21  +  C  sin  21  ]  , 

vz; 


z  z 

where  M  =  J  p(z')dz',and  I  =  J  KQ(z')dz'.  The  coefficients  A,  B,  C  can  be 
o  o 

found  by  using  the  initial  conditions  for  a  matched  beam,  a(z=0)  =  a  , 

1  »  » 

a  (z=0)  =0,  a  (z=0)  =  0  .  The  matched  beam  radius  aQ  is  related  to  the 

H  2  2  2 

initial  transverse  emittance  a  =  e  (0)/(K_(0)Y  (0))  .  Using  the  initial 

o  n  B 

conditions,  we  arrive  at  the  following  expression  for  the  root-mean-square 
beam  radius, 


a  =  a  e 
o 


-M/2 


V0)  il  /2 

kb(z) 


[  1 


u(0)  +  Kg(0 ) /Kg(0 ) 

2Kb(0T 


sin 


21  ] 


1/2 


(13) 


Equation  (13)  shows  that  the  beam  radius  does  not  remain  constant  even  when 
the  beam  is  matched  at  injection.  In  addition  to  the  exponential  decay  from 
the  radiation  damping,  the  beam  envelope  developes  an  induced  betatron 
oscillation.  However,  the  normalized  emittance  is  just  an  exponential  decay 
given  by  Eq.  (11a). 

We  may  gain  some  insight  into  the  general  effect  of  radiation  damping  on 
the  transverse  emittance  by  studying  Eq.  (12)  in  the  case  of  untapered  wiggler 
field.  We  shall  first  consider  the  case  where  Y’  =  0.  This  could  be  the 
situation  when  the  acceleration  mechanism  is  saturated  by  the  radiation 
damping  and  the  beam  energy  is  constant.  The  evolution  of  the  beam  radius  is 
then  given  by  the  appropriate  limit  of  Eq.  (13).  Since  there  is  no  tapering 


of  the  wiggler,  the  solution  is  exact  and  given  by 


a 


/  C.  f  ,  .  V  .  ~ I  1C. 

V  [  1  2KT  Sin  2KBZ  ] 


V 


The  beam  radius  again  exponentially  decays  with  an  induced  betatron 

oscillation.  Since  Y  is  constant,  the  damping  rate  v  is  constant,  and  the 

normalized  emittance  e  is  given  by  e  (z)  =  t  (0)  exp(-vz). 

n  n  n 

Next,  we  consider  the  situation  when  an  accelerated  beam  is  cooled  by 
passing  it  through  an  untapered  external  wiggler  field.  Since  the  beam 
decelerates  due  to  the  synchrotron  radiation  damping,  we  have  Y'/Y  =  -v. 

This  gives  p  =  0  and  since  K_  =  a  k  /(/2Y),  the  betatron  wave  number  KD  is  a 

D  W  W  D 

function  of  z.  The  spatial  dependence  of  Y  can  be  evaluated  using  Y'/Y  =  -v, 

2  2  2  2 

and  Eq.  (13)  reduces  to  a  =  a  ( 1  +  v  z  ),  where  v  =  t0a  k  Y  c.  Although  the 

oo  oRwwo 

p 

beam  radius  remains  constant  up  to  order  of  0(z  ),  the  normalized  beam 

emittance  decreases  algebraically,  e  =  e  (0)/(1+v  z). 

n  n  o 

The  relevance  of  the  above  analysis  depends  on  the  magnitude  of  the 

p 

damping  rate  vq.  For  the  following  set  of  accelerator  parameters, 

El  =  1.5x10^  V/cm.,  Bw  =  50  kG. ,  A  «  1  m,  it  is  estimated  that  the  e-fold 

5 

length,  1/vQ,  could  be  as  short  as  <  600  m  for  Yq  =  10  .  Therefore,  our 
results  show  that  one  can  improve,  by  induced  synchrotron  radiation,  the 
quality  of  an  electron  beam  by  passing  it  through  an  external  wiggler  field. 

I V .  Quant urn  Excitation 

An  estimate  for  the  minimum  transverse  normalized  beam  emittance  due  to 

quantum  excitation  in  an  IFEL  accelerator  can  be  obtained  from  the  following 

qualitative  treatment.  Similar  arguments  can  be  made  for  electron  beams  in 
17  1 

storage  rings.  ''  The  normalized  transverse  velocity  and  radial 

displacement  of  an  electron  in  a  wiggler  field  are  given  by  Bw=aw/Y,  and 

r  =a  \  /(2itY).  For  a  fluctuation  5E  in  the  energy  of  the  electron,  the 
w  w  w 

corresponding  fluctuations  in  r  and  8  are  6r  =n,5E/E,  and  68  =£6E/E,  where 

w  w  w  w 


MrM*  W 


10 


n=a  X  /(2itY)  and  £= a  /Y.  The  increase  in  normalized  emittance  due  to  such 
w  w  w 

fluctuations  is1^’1^  Ae  =  Y[KD<6r2>  +  <6B2>/KD],  which  for  a  weakly  focusing 

n  d  w  w  b 

channel,  KR<<k1.,  can  be  approximated  by  Ae  *  Y<6B2>/KD  -  (Y£2/Kd)<6E2>/E2 
D  w  n  W  D  D 

2 

Due  to  the  discrete  nature  of  the  synchrotron  radiation,  <6E  >  is  given 
2 

by  N(rtiu  )  where  N=Pz/(ctfio  )  is  the  number  of  photons  emitted  in  a  distance  z 
c  c 

P  is  the  synchrotron  radiation  power,  and  is  the  energy  associated  with  a 
quantum  of  synchrotron  radiation.  We  can  therefore  obtain  the  rate  of  change 
of  en  due  to  quantum  excitation, 


/V  Y£f 

d*  «-E-'  kb  cE2  ‘ 

However,  with  radiation  damping,  the  total  change  in  is  given  by 


,  n,  .  /  n » 

(~ — )  »  -  ve  +  (- — )_  _ 
dz  n  dz  Q.E. 


The  normalized  emittance,  e  ,  reaches  a  minimum,  de  /dz=0,  when  the  two 

n’  n 

effects  are  balanced.  This  gives  e  =Y £2hw  /(KDE)  for  the  minimum  normalized 

n  c  d 

emittance,  where  we  have  used  vc=P/E.  For  synchrotron  radiation, 

r1oj  =3rtcY^/(2p)  where  p  =  Y/(a  k  )  is  the  radius  of  curvature  of  the  electron 
c  w  w 

orbit  in  the  wiggler.  The  minimum  transverse  normalized  beam  emittance  is 


then  approximately  given  by 


(e  )  .  =  3fiaJk  /(2m  cKn)  . 

n  min  w  w  o  B 


In  the  case  of  weak  focusing  due  to  wiggler  transverse  gradients, 

K  =  a  k  /(/2Y),  and  the  minimum  normalized  emittance  is 
B  w  w 


(e  )  . 
n  min 


3tf Ya  / (/2m  c)  . 
w  o 


Using  the  accelerator  parameters  at  the  end  of  section  III,  Eq.  (15) 
gives  the  value  of  the  minimum  normalized  emittance  to  be  -1.8  cm-rad. 

Such  a  large  value  of  the  minimum  emittance  indicates  the  inadequacy  of  the 
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weak  focusing  from  the  wiggler  transverse  gradients.  Strong  focusing  from. 


for  example,  a  rotating  quadrupole  field  produced  by  a  pair  of  (or  four) 

prv  p  i 

helical  current  windings  ’  may  be  required.  The  betatron  wavenumber  for 

??  2  2 

such  a  focusing  mechanism  is  given  by  KD  =  lei ( 9B/3r )/Ym  c  ,  where  9B/3r  is 

D  O 

the  magnetic  field  gradient  of  the  quadrupole  field  on  axis.  For 
3B/3r-250  G/cm,  a, .=600,  \  =10m,  and  Y=4x10  ’  Eq.(l4)  gives  a  minimum 

W  w 

normalized  emittance  of  e  -0.13  cm-rad.  Another  possible  strong  focusing 

force  could  be  the  electrostatic  radial  electric  field  of  an  ion  column.  Such 

a  column  could  be  created  by  the  ionization  of  the  residual  gas  by  a  low 

energy,  high  current  electron  beam  pulse  preceeding  the  main  accelerating  beam 

pul.'e.2^  The  betatron  wavenumber  for  such  a  focusing  mechanism  can  be 

2  2  2 

easily  shown  to  be  K  =  m  . (m./m  )/(2Yc  ),  where  w  .  is  the  ion  plasma 

d  pi  i  o  pi 

frequency  and  (m^/m0)  is  the  mass  ratio  between  the  ions  and  the  electrons. 

For  n1  =  1012/cm^,  aw=600,  Xw=10m,  and  Y=4x105,  Eq.  (14)  gives  a  minimum 

normalized  emittance  of  e  -0.04  cm-rad.  An  additional  benefit  of  having  ion 

n 

focusing  in  the  IFEL  accelerator  is  that  the  radial  plasma  electron  density 
profile  in  an  ion  column  can  also  be  a  focusing  medium  for  the  laser  beam. 


3 


where  R1  =  /2  |  e ]  ELsim|>R/ (m0c2/k)  ,  R2  =  | e | ^k/( Sm^c 1 0 )  ,  is  the  resonance 
phase,  E^  the  laser  electric  field  strength,  and  k  the  laser  wave  number.  The 
second  condition  together  with  the  pendulum  equation,  Eq.  (6),  provide  the 
spatial  dependences  of  kw  and  B  , 


2k  4E.  k^ 

3  k;  -  nr  b:  +  —  sin  *r  r 

w  w 


2v'*2mc* 


:  |/k 


/r 

w 


=  0. 


(17) 


Equation  (17)  shows  that  the  required  tapering  of  the  wiggler  field  may 
be  obtained  by  prescribing  and  a  relationship  between  kw  and  Bw  in 
Eq.  (17).  As  an  example,  we  assume  the  tapering  of  the  wiggler  field  to  be 

p 

that  of  a  maximum  rate  IFEL  accelerator.  For  such  a  case  the  wiggler 
strength  and  the  wiggler  period  are  related  by  the  following  power  law, 


(R1/6R2) 


\/k  ,  7/8 


Equation  (17)  may  then  be  solved  to  give 


where 


Bw  =  bw«»[i  +  V]-7/9* 


kw  =  kw(0)[1  +  V] 


-8/9 


(18a) 

(18b) 


R  -  —  R  (R  /6R  )3/7B  (0)9/7  . 

“  |e|A  2  '  2 

Evaluating  Eq.  (11)  and  (16a)  with  (16c)  and  (I8a,b)  gives  the  normalized 
transverse  emittance  and  the  resonant  energy  of  the  beam  as  functions  of  the 


propagation  distances. 


For  our  example,  we  will  consider  the  following  set  of  accelerator 
parameters2:  E.  =  1.5x10^  V/cm,  B,,(0)  =  50  kG,  A  (0)  =  100  cm,  X  =  10.6  pm 
and  a  resonance  phase  of  sin  =  0.6.  The  initial  conditions  are  for  a 

matched  beam  with  a  radius  of  1  mm  ,  a  normalized  emittance  of 
eQ*.205  cm-rad,  and  the  required  beam  injection  energy  is  -  52  GeV .  The  beam 
is  allowed  to  propagate  for  1  km  without  depleting  the  laser  radiation.  We 
repeated  the  calculation  by  assuming  there  is  no  radiation  damping  but  with 
the  same  power  law  tapering  of  the  wiggler  field. 

The  results  are  represented  in  Figs.  (1),  (2),  and  (3).  The  open  squares 
denote  the  presence  of  '-adiation  damping,  while  open  circles  denote  its 
absence.  From  Fig.  (1),  we  can  see  that  the  final  energy  is  not  significantly 
reduced  by  the  radiation  damping.  Figure  (;? )  snows  the  exponential  decay  of 
the  normalized  emittance.  At  the  ~nd  jf  the  t ne-<i lometer  accelerator,  the 
normalized  emittance  is  reduced  t  ei-rad  which  is  very  close  to  the 

minimum  normalized  emittance  of  -0.  ••  im-rad  at  th  it  point  if  ion  column 
focusing  is  assumed  in  the  iccelerat or .  In  Fig.  (3',  the  appropriate  tapering 
of  kw  and  Bw  for  the  two  eases  are  shown. 

Conclusion 

We  have  studied  the  evolution  of  transverse  emittance  and  the  beam  radius 
due  to  the  radiation  damping  effect  in  an  1FEL  accelerator.  We  derived  the 
beam  envelope  equation,  Eq.  (10),  which  includes  the  effects  of  radiation 
damping,  and  have  demonstrated  that  the  normalized  transverse  emittance 
decreases  exponentially  with  a  damping  rate  given  by  the  radiation  damping 
coefficient  v  until  it  reaches  a  minimum  value  due  to  quantum  excitation 
The  beam  envelope  equation  was  solved  analytically  for  a  slowly-varying 
wiggler  field.  We  have  derived  an  expression  for  the  minimum  normalized 
emittance  in  the  IFEL  accelerator  and  showed  that  strong  focusing  is  essential 


in  reducing  this  minimum  emittance  due  to  quantum  excitation.  We  have  shown 
that  radiation  damping  can  play  an  important  role  in  improving  beam  quality 
without  a  significant  sacrifice  in  beam  energy. 
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